Gradient corrections for semiclassical 
theories of atoms in 
strong magnetic fields 

Christian Hainzl 

Institut fiir Theoretische Physik, Universitat Wien 
Boltzmanngasse 5, A-1090 Vienna, Austria 

E-Mail: hainzlSthp . univie . ac . at 

November 29, 2000 

Abstract 

This paper is divided into two parts. In the first one the von Weizsacker 
term is introduced to the Magnetic TF theory and the resulting MTFW 
functional is mathematically analyzed. In particular, it is shown that the 
von Weizsacker term produces the Scott correction up to magnetic fields 
of order B <^ Z^, in accordance with a result of V. Ivrii on the quantum 
mechanical ground state energy. 

The second part is dedicated to gradient corrections for semiclassical 
theories of atoms restricted to electrons in the lowest Landau band. We 
consider modifications of the Thomas-Fermi theory for strong magnetic 
fields (STF), i.e. for B <^ Z^. The main modification consists in replacing 
the integration over the variables perpendicular to the field by an expan- 
sion in angular momentum eigenfunctions in the lowest Landau band. This 
leads to a functional (DSTF) depending on a sequence of one-dimensional 
densities. For a one-dimensional Fermi gas the analogue of a Weizsacker 
correction has a negative sign and we discuss the corresponding modifica- 
tion of the DSTF functional. 
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1 Introduction 

In this paper we study gradient correction terms for semiclassical theories describ- 
ing the ground state energies of heavy atoms in strong homogeneous magnetic 
fields. Such an atom, with N electrons of charge — e and mass rUe and nuclear 
charge Ze is described by the nonrelativistic Pauli Hamiltonian 

i<i<A^ I il J l<i<j<N ' * ^' 

acting on the Hilbertspace L^(R^, C^) of electron wave functions. The 

units are chosen such that h = 2me = e = 1. The magnetic field is B = (0, 0, B), 
with vector potential A = ^B{—X2,Xi,0), where B is the magnitude of the field 

2 3 

in units of Bq = "^^3 ^ = 2.35-lO^Gauss, the field strength for which the cyclotron 
radius Ib = {he/ (ei?))^/^ is equal to the Bohr radius oq = / (meC^). The ground 
state energy is 

E^{N, Z, B) = inf {(^, Hn^Ij) : V e domain Hn, (t/-, V^) = 1}. (1.2) 



In | |LS Y2 | Lieb, Solovej and Yngvason approximated ( |1.2| ) by means of the 



MTF (Magnetic Thomas-Fermi) functional 

S'^'^^lp] = J rBip)-J Vp + Dip,p), (1.3) 

where V{x) = Z/\x\ and D{p,p) = |(p, * p). The magnetic energy density 
tb is, by definition, the Legendre transform of the pressure Pb, i.e. 

TB{t) = snp[tw - PBiw)], (1.4) 

with 

R 00 

Psiw) = T^(«;'/' - 2 V \2tB - wfj'). (1.5) 

The corresponding energy 

E^^^{N,Z,B) = inf{^^T^[p]|p > 0,p G D^^"" , j p < N} (1.6) 

was proved by Lieb, Solovej and Yngvason to be asymptotically exact, as shown 
in the following Theorem: 



1.1. THEOREM, f ^n^m Theorem 5.1) If Z ^ 00 with N/Z fixed andB/Z^ 
0, then 

E'^iN, Z, B) /E^^^{N, Z, B) 1. 
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In the hmit B oo the function tb is the kinetic energy describing particles 
confined to the lowest Landau band, i.e., 

which results in the STF (Strong Thomas-Fermi) functional 

S'''^[p] = ^,Jp'-Jvp + D{p,p). (1.7) 

The corresponding energy E^"^^ is quantum mechanically exact in the limit Z — >■ 
oo for Z'^^^ <ti B <^ Z^ (| [LS Y2| Proposition 4.16), which emphasizes the fact 
( IILSYlj Theorem 1.2) that for B > Z^/^ ^j^g electrons are to leading order 
confined to the lowest Landau band. 



1.1 Corrections to the leading order of the full Hamilto- 
nian ( |1 . 1|) 

The best result to date concerning corrections to the leading order of ( |1.2| ) is 
presented by Victor Ivrii in m. Theorem 0.2: 

1.2. THEOREM, (gj Theorem 0.2) Let B < Z^ and N ^ Z, then 

|E«(Ar, Z, B) - E^T^(Ar, Z, B) -^Z^\<Ri + i?2, (1-8) 

with 

Ri = CZ^/\N + Bf/^ and R2 = CZ^I^B^'^. (1.9) 

Recall the order of the energy. 

We make a few comments concerning Ivrii's proof. Let 

i^A = [(-^V + A(x)) ■ (t]2 (1.10) 

denote the free Pauli-Hamiltonian and 

0^T^ = z|xri-p^T%ixri (1.11) 

the self consistent magnetic TF potential. The main part of the estimate ( |1.8|) is 
given by the difference between 

Tr[i/A-0^^^ + ^]-, (1.12) 

the sum of all negative eigenvalues of the operator Hp^ — (^"^^ + /i, and its 
semiclassical approximation 

jpsir^'-p). (1.13) 
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For those who are famihar with microlocal analysis we should remark that Ivrii 
does not really consider 0^^^ but a smooth mollification, which we also denote 
with 0^"""^ for simplicity. In order to derive accurate estimates of 



Tr[/JA - 



iMTF 



iMTF 



;i.i4) 



Ivrii essentially divides the domain into two main zones, for B > Z^^^, namely 

Xi = {x| < |x| < B/Z} and X2 = {x| B/Z < |x| <rs = Z^^^B'^^^}. (1.15) 

A corresponding partition of unity is given by two function ipi and !f2, with 
^i + f2 = 1 on X1UX2 and ipi essentially supported in Xi- Using scaling arguments 
and semiclassical spectral asymptotics Ivrii treats each zone separately. In the 
inner zone xi? where all Landau levels are taken into account and the MTF 
potential is very similar to the usual TP potential, he gets 



Xi ■ 



Tr(^i[i/A 



iMTF 



iMTF 



12) - -Z^ 
4 



<i?i. (1.16) 



We see that in xi the Scott correction is recovered. Moreover, we should note 
that the machinery of semiclassical spectral asymptotics can only be applied to 
Xi under the condition Z/B ^ 1/Z which means that ( |1.16|) is only valid for 
B < Z"^^^, with arbitrary 6 > 0. For B > Z^ a. semiclassical approximation is no 
longer possible and the terms ( |1 . 1 2|) and ( |1.13|) have to be estimated separately. 
Since R2 overcomes Z^ for B > Z^/^, we should point out that the Scott correction 
in ( |1.8| ) only provides the next to leading order for B <C Z^/^, but in the domain 
Xi it nevertheless makes sense up to S <C Z^ according to ( [1.161 ). 

In the outer zone X2 only the lowest Landau band is occupied, which implies 
that in this region the MTF energy is represented by the STF energy correspond- 
ing to the functional (|1.7|) . In X2 Ivrii derives the estimate 



X2 



-MTF 



f2PE 



iMTF 



< R2 



;i.i7) 



where the main contribution of ( |1.17|) really stems from the edge of the STF atom 
rs ~ 

For low magnetic fields {B < Z) V. Ivrii even improves ( |1.8|) and recovers 
Dirac and Schwinger corrections as well. 



1.3. THEOREM, (gj Theorem 0.3) If B < Z then 
EQ^(iV, Z, B) - E^TF(^^ B) - -^Z^ + CDS j {p^'^Y'^ 

holds with some (5 > and an appropriate parameter Cds- 



< cn^/'\{i + b)/nY 

(1.18) 
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1.1.1 The von Weizsacker term introduced to MTF 

The von Weizsacker correction term was successfully introduced to the TF theory 
in the sense that it reproduces the Scott correction (rigorously proven in ||SW| , 0), 

E™ = ET^ + 0(Z2) + o(Z2). (1.19) 

In addition to the correction, the TFW theory remedies some defects of 
the TF theory: The corresponding TFW density is finite at the nuclei, bind- 
ing of atoms occurs and negative ions are stable, furthermore the density has 
exponential fall off at infinity, at least for neutral atoms and molecules. 

In the limit B ^ the function tb is the kinetic energy density in zero 
magnetic field, i.e. 

ro(t) = -(37r2)2/3t5/3_ 
5 

Since for small values of B the introduction of the von Weizsacker term to the 
MTF functional is justified, we will further check up to which values of B this 
definition makes sense. We get the functional 

£:™[p] =A J |Vpir + I Tsip)- lvp + D{p,p), (1.20) 

with a suitably chosen parameter A. The functional ( |1.20| ) can be treated anal- 
ogously to the usual TFW functional in ||BBL||. So we will just sketch the proofs 



of the main propositions. 

It turns out that for B <^ the von Weizsacker term still produces the Scott 
correction, but makes no longer sense for higher magnetic fields. We will derive 
the following Theorem: 

1.4. THEOREM. For all B, Z and N/Z fixed 

I^MTFW _ ^MTF _ Q^^2)| < ^'5^/5^3/5 _^ ^(^2) ^^ ^l) 



Remark: The estimate ( |1.21D is clearly useful if i? < Z"^^^. Theorem |1.4| will 
be proved in Section 2.2. 

The Scott correction, just like in TFW theory, comes from distances of order 
1/Z near the nucleus, whereas the bound CB^/^Z^^^ comes from the edge of the 
MTF atom and dominates the Scott correction for B > Z"^!^. Moreover, ( 1.21| ) 



is in accordance with Theorem |1.2| , which justifies a posteriori the introduction 
of the von Weizsacker term to the MTF functional. 

1.2 Physics in the lowest Landau band 

The quantum mechanical ground state energy of particles confined to the lowest 
Landau band is given by 

Elf(iV, Z, B) = inf (^P, U^Hr^U^^P). (1.22) 

IIV'l|2=l 
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where Ho represents the projector on the lowest Landau band, given by the kernel 
Ilo{x,x') = ^exp x ■ B - -x'_^fB^ 6{x3 - X3)P|, (1.23) 

where Pj denotes the projection onto the spin down component, and 11^ denotes 
the A^'th tensorial power of Ho. The leading order of -E^nf ' SiS B, Z ^ oo with 
B <^ Z^, is given by E^"^^, the ground state energy of the functional ( |1.7|) . In a 



companion work we show, what is expected by Ivrii's Theorem |L2 



I^Q^^ _ ^STF| < ^54/5^3/5_ (1 24) 



As we have argued above the main contribution to the estimate ( p..24| ) comes 



from the edge of the STF atom, rs ~ Z^^^B Recall the order of E^^^ , 
ESTF[iV,Z,P] = Z^{B/Z^f/^E^^^[N/Z, 1,1], so that the estimate is only of 
interest for B <^ Z^. 

As a better approximation to -E^^f, valid also for B > Z^, Lieb, Solovej and 
Yngvason suggested a density matrix functional defined as 

^°^[r]= / Tr^2(K)[-9|r,.Jdxx-Z / \x\-'pr{x)+D{pr,pr). 
Jr2 J 

Its variable is an operator valued function 

T : x±^ r^.^, 

where F^^ is an integral operator on L^(R), given by a kernel rx^{x3,y3) and 
satisfying 

< < (P/27r)/ (1.25) 
as an operator on L^(R). The energy 

E^^{N, Z, B) = inf{^°"[r]| r satisfies ( [i:25|) and j pv < N} 

turns out to be asymptotically exact for magnetic fields in the following precise 
sense: 



1.5. THEOREM. (^^SY^ Theorem 5.1 and 7.1 and equations 1.3 and 8.5) 
For some constants, C\ and C^, we have 

Ru > (AT, Z, B) - E^^{N, Z, B) > ~Rl, (1.26) 

with 

Rl = CAmin{Zi^/i^p2/5^^8/3^^ ^ (ln(P/Z3))2]} (1.27) 

and 

Ru = C'^mm{Z^/''B^'\Z'"'^[l + \n{Z) + (ln(P/Z3))2]5/6}. 
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We remark that the STF energy is the natural semiclassical approximation of 
the DM energy. More precisely, the DM energy can be written as 



B 
2^ 



d^^TT,.^^)[-dl - - + /.°^iV - p°^), (1.28) 



whereas the STF energy is given by the corresponding semiclassical expression 

^ y rfx^ y j dpdz\p' - - /^'^^]- + /^'^^^ - ^(p'^^ p'^^)- (1-29) 

With the decomposition L^(R^ dx; C^) = L'^iR'^, dx±) (g> L'^(R, dz) ® the pro- 
jector Ho can be written as 

no = A, (1.30) 

m>0 

where ipm denotes the function in the lowest Landau band with angular momen- 
tum —m < 0, i.e., using polar coordinates {r,ip), 



B 1 fBr 



2\ ™/2 



'''-> = ^/^;^l^l e-'-e--/'. (1,31) 



Using this and H^^rn = 0, we can write 



m>0 



Based on this decomposition the author and R. Seiringer introduced in |PS | a nat- 



ural modification of the DM functional called discrete density matrix functional 
(DDM) 

^bT\^] = E (Tr[-^'r^] - ^ / VUz)Pm{z)dz^ + D{p,p), (1.33) 



mSNo 

where 



-D(p,p) = ^E / Vm,n{z - Z)pm{z)pn{z)dzdz , (1.34) 

and the potentials Vm and Vm,n are given by 

Kn(2;) = j j^|</>m(x_L)|^(ix_L, 

Vm,n(2--2) = / , 71 d^±d^^. (1.35) 
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Here F is a sequence of fermionic density matrices acting on L^(M, dz), 

r = (rm)m6No5 (1.36) 

with corresponding densities p = {pm)m, Pm{z) = Tjn{z,z). Note that 
depends on B via the potentials Vm and Kn,„. The corresponding energy is given 
by 



E™(Ar, Z, B) = inf <^ S^Ti^] 



5^Tr[r„] <iV I. (1.37) 



It is shown in [HS| that E^^^ correctly reproduces the confined ground state 
energy -E^^f ^P^-^'t from errors due to the indirect part of the Coulomb interaction 
energy: 



1.6. THEOREM. f HH^J Theorem 1.2) For some constant c\ depending only on 
X = N/Z 

> (iV, Z, B) - E^^^{N, Z, B) > -Rl, (1.38) 

with 

Rl = cxmm{Z^'/^^B^/^ Z^/%1 + [lniB/Z^)f)} . (1.39) 

Since the functional ( |1.33| ) can also be seen as a reduced Hartree-Fock func- 
tional, in the sense of |S|, it does not surprise that the upper bound in (|1.38| ) is 



an improvement to ( |1.26| ), the relation between E^^^ and -E'^nf In addition to 



better estimates, the DDM theory remedies the defect of the DM theory having 
a sharply cut ground state density supported in the set {x| |x| < ^/2N/ B}, for 
the respective three dimensional DDM density, 

P°°^(x) = $:P™°''W'^™(-^)' (1-40) 

m 

has exponential fall off at infinity. Furthermore, since the DDM energy describes 
-^ranf correctly apart from errors due to the indirect part of the Coulomb interac- 
tion energy, the DDM energy could give rise to even recover the exchange term, 
by means of an improved lower bound on the two body Coulomb repulsion for 
particles in the lowest Landau band. For the exchange energy is anticipated to 
be of order ln(5/Z^/3) ^7/5^1/5 B <^Z^, which in fHSf by the author and R. 
Seiringer is conjectured to be given by the term 

cln(i?/Z^/3)5^ ApP°^)^ (1.41) 
i ^ 

This would lead to the relation 

^Q^^ = i^DDM _ cin(i?/Z^/3) ^ j (ppDM)2 ^ ^ (\n{B / Z^/'') Z'/' B^/') , (1.42) 
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with c appropriately chosen. 

We have stated above that the STF functional is the natural semiclassical 
approximation of the DM functional. Hence, we can ask for the natural semiclas- 
sical approximation of the DDM functional, of which the answer is given by the 
so called DSTF functional 



^ ^ I pUz)-Z I V^{z)p^{z)dz]+D{p,p), (1.43) 

meNo 

where p is a sequence of one-dimensional densities, p = (pm)meNo5 = ^r^/S, and 
the respective DSTF energy is defined as 



^DSTF^^^ Z,B) =inf <^£°STF[^] 



m J 



;i.44) 



In Section 3.2.2 we will argue that the DSTF functional is even the natural 
semiclassical approximation of the ground state energy -E^^f itself. 



1.3 Gradient corrections for semiclassical lowest Landau 
band theories 

1.3.1 The Tomishima-Shinjo correction term 

For higher magnetic fields, where only the lowest Landau band has to be taken 
into account, Tomishima and Shinjo |P?S|| obtained the gradient correction term 



I.e. 



^^'[p] = |r / Pi^^p)' - 1 1 (^\\p'"f + (1-46) 

by perturbation expansion of the canonical density matrix. In 1995 the authors of 



[MZP|| recovered the TS theory within the framework of current density functional 
theory. Since ( |1.46| ) has a negative gradient correction along the magnetic field 
the TS functional is no longer bounded from below. Hence the corresponding 
energy cannot, as usual, be defined by minimizing over a suitably domain of 
definition, but only through the solutions of the corresponding Euler-Lagrange 
equation under the restriction f p = N, i.e 



^(V„P)^ - A,p 



V - p* ^ - p(N). 

\x 

(1.47) 
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A direct attack on this complicated equation does not look promising, but a rough 
estimate of the corrections to STF can be obtained by inserting the density 



into (|1.46| ). The negative gradient term gives a correction -0{B^/^Z^/^) commg 



from the edge of the STF atom. From ( p..24| ) we know that 

I^Q^^ _ ^STF| < CB'/'Z^/'', (1.49) 

where the main contribution also stems from ~ Z^^^B~'^^^, the edge of the 
STF atom. 

On the other hand the positive gradient correction orthogonal to the magnetic 
field in ( |1.46| ) produces a correction O^B^^'^Z^^'^) at a distance of order i?"^/^ from 
the nucleus. This part of the correction can also be obtained from an isotropic 
Tomishima functional, defined as 

S'^[p]=S^^^[p] + ^-^ j p{Vpf. (1.50) 



The functional ( |1.50| ) has all the good properties of the usual TF theories, such 



as convexity and boundedness from below. The study of this functional, which 



we do in detail in Section 3.1, should help us to get a deeper understanding of 



the nature of the positive correction term in (|1.46| ). For the ground state energy 



of ( |1.50| ) we will derive the following theorem. 



1.7. THEOREM. For all B, Z and N/Z fixed 

E'^{N, Z, B) - ESTF(iv, Z, B) = OiB^'^Z""^) + o{B^'^Z'"^). (1.51) 



Furthermore we will argue in Section ^]T] that the p{W pY term remedies the 
defect of the STF theory that the full Coulomb potential is used although the 
particles in the lowest Landau band do not see the full singularity, since they are 
smeared over a region of radius B^^^"^. In contrast to TFW theory, where the 
maximal number A^^ of electrons that can be bound is strictly larger than Z, it 
will as a slight surprise turn out that in IT theory = Z, just like in the STF 
theory itself. Also, the radius of atoms in IT theory is finite, as in STF theory. 
These features confirm that the p(Vp)^ term essentially only effects the density 
close to the nucleus. 



1.3.2 Gradient correction for the discrete STF theory 

As discussed above, the gradient term ~ p(V_Lp)^ in ( |1.46D produces essentially 
a smearing of the Coulomb singularities over a distance of the radius B"^/^. The 
same effect was obtained by replacing STF by DSTF. 
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It appears thus natural to look for a negative gradient term that has an 



analogue effect in DSTF as the negative gradient term in (|1.46|) has in STF 
theory, i.e. provides corrections at the edge of the atom. 

The DSTF theory is effectively a theory of coupled one-dimensional problems. 
Analogous arguments as lead to the von Weizsacker term for a three dimensional 
Fermi gas give for a one-dimensional Fermi gas a gradient correction — ^(Vp^/^)^ , 



cf. [ph| . Hence we suggest the definition of a discrete von Weizsacker iunctional: 



£^^[p] =J2(-ll \d.^/pU^)\' + ^ I pUz)-Z j V„.{z)p„,{z)\ +D{p,p) 

meNo 

(1.52) 

By denoting ^/p^ = ipm we arrive, under the restriction I V^m ~ 
corresponding TF equation 

(l/3)9,Vn(^) + 3«V^(z) = [ipniz) - /i(iV)]^„(z) Vn e No, (1.53) 

with 

= ZVniz) -Y,j '4^l{z')V^,n{z - z')dz' . 



These coupled equations are probably somewhat easier to deal with than ( |1.47|) . 



If we reduce ( [1.52| ) to the angular momentum channel m = and drop the 
Coulomb repulsion term, we get the one- dimensional functional 

^'^^[p] = -\j \d.^)? + j pHz) -Z j V,{z)p{z). (1.54) 

This simplified functional will be studied in Section 3.3.2. In particular, we shall 
show that the negative gradient term reproduces the right QM correction to the 
energy without the gradient term. 

2 The magnetic TFW theory 

2.1 Mathematical analysis of the MTFW functional 

In this section we are going to mathematically analyze the MTFW functional 

£™[p]=A j |Vpip + j tM- Jvp + D{p,p). (2.1) 

Since the mathematical propositions do not depend on the parameter A, we let 
A be 1 in this section. The most important features of TB{t) which will be used 
in our calculations are (compare |[LSY2|| Lemma 4.1): 



3 

T'Jt) < Kit^/^ and TBit) < (2.2) 

5 
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with Ki = (47r2)2/3. 

Since the MTFW functional does not differ very much from the functional 
BBL||, where the authors used a kinetic energy density r(p) = {l/p)ff, our 



m 



procedure in analyzing (|1.2CI|) will be in analogy to their work. We are thus 
concerned with the minimizing problem 

Min{£^™[p]|p > 0,p e n Lff , Vp^ G and j p = iV}, (2.3) 

where is a positive constant, which physically is the total charge number. Our 
main result is the following: 

2.1. THEOREM. There is a critical number < Nc < oo, so that 



1. if N < Nc ( \2.!^ ) has a unique minimizer, 



2. if N > Nc ( \2.^ ) has no minimizer, 

3. Nc > Z. 



Similar to |PBL|| we first examine the problem 

Min{£^™[p]|pGD} (2.4) 

with 

D = {p|p > 0, rsip) <oo,pe L\ Vp^/2 G L\ D{p, p) < oo} (2.5) 

and proof the existence of a unique minimizer po. Since D contains the domain of 
( |2.3| ) we have to show po G L^(]R^) in order to guarantee Nc < oo. Furthermore 
we will derive Nc > Z, which shows that this theory allows negative ions. The 
proofs will be based on the Euler-Lagrange equation for ip = ^/p^. 
First we consider some basic properties of (p.l|) . 



2.2. LEMMA. For D defined in we have 



D C {p|p > 0, p G L3 n L\il ^ jj^p^ p) < oo} = D. (2.6) 

Proof. According to ||LSY2|| (4.19) one gets for all 9.(1^ 

[ p{xfl^dx < — [ TB{p{x))dx + C VoK^a) < oo 

with Q = QiUQ2- □ 

2.3. PROPOSITION. The absolute minimum of £^^^^[p] is achieved for a 
unique po E D. 
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Proof, (cf. [BBL] Lemmas 2,3,4 and 5.) Recall that, by definition ( p..4| ), TB{t) is 
strictly convex, hence £^^'^^^[p] is strictly convex. 

Let pn be a minimizing sequence. There exists a constant C such that 

llPnIla <C, [ TB{pn) < C, \\V p]!% < C,D{p^,pn) < C. 



By the Banach-Alaoglu theorem we can extract a subsequence, still denoted as 
Pn, with 

Pn Po weakly in L^, (2.7) 
Vpl/^ Vp]J^ weakly in L^. (2.8) 

Since by use of Holder's inequality ||py^||_H-i(n) < C{VL) and H^{VL) is relatively 
compact in I/^(n), if is a bounded smooth domain, pH"^ has a subsequence 
converging in LP'{Vt). Using Cantor's diagonal trick on a sequence of increasing 
r2's we arrive at 

pT^pT a. e.. (2.9) 

By Fatou's Lemma we get 

liminf j tb{p„) > j rB{po) and liminf L'(p„, p„) > L'(po, po). 
Since U norms are weakly lower semicontinuous, 

liminf j iVp^T > / iVp^ p. 
Moreover, one can show, in analogy to Proposition |3.3|, that 



so we altogether arrive at 

liminf £^™[p„] > ^^™[po]. (2.10) 

The uniqueness follows from the strict convexity of £^'^"'"^^[p]. □ 

For the minimizing po we now can derive an Euler-Lagrange equation. Denote 

2.4. PROPOSITION. The minimizing ip"^ = po satisfies 

-Atp + T's{4j')^ = ^^, (2.11) 
in the sense of distributions, with (f = V — ip"^ * . 
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Proof, (cf. [BBL] Lemma 6.) Note that ip"^ & D implies ipip , t'q{iIj'^)iIj G Ll^^^ 
which gives (|2.11|) a meaning in the sense of distributions. Consider the set 



D = {C\CeL^n , VC e and D{C\ e) < oo}. 
If C G ^ then p = C e D a.nd 



(2.12) 



We find < 0(C) for all C e ^- Let G . Using the fact that |0(?/^ + 
tvi)\t=Q = 0, we easily arrive at 



ipArj + / Tq{iIj )il)ri = / (fipt]. 



(2.13) 



□ 



Starting from Equation ( |2.1lD we can now step by step gain several properties 
for ijj. 

2.5. LEMMA, ip is continuous on M^, more precisely G C[°'" for all a < 1. 
Proof (cf. [BBL] Lemma 7.) Since (|2Tl| ) yields -Aip < ip^j, with ipip G Lf'J, 



one gets G LJ^^. Again using (|2.11|) the proposition follows by means of standard 
elliptic regularity theory. □ 



2.6. PROPOSITION. ^ G L 



Proof, (cf. [[BBLf Lemma 8.) Assume, by contradiction, J ip"^ = oo. Then we can 
choose an r, such that 



^^(x) > Z + 25, 



|x|<r 



for some 5 > 0. Therefore 



ip'^*\x\-^> / V^^(x)(|x| + li/D^^c/y > (Z + 25)/(|x| +r). 



which gives us 



if{x) = V{x) -ip"^ *\x\ ^ < 



Z + 25 



\x\ — r \x\ + r 
with |a;| > r. Thus there exists an ri > r, such that for |x| > ri 



(2.14) 
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From ( p.ll| ) we get 

- AV' + 5|x|-V < 0, (2.15) 
for |x| > ri. Now we choose a comparison density 

V5(x) = Me-2(W/^, 

which satisfies 

- Ai> + S\x\-^i^>Q. (2.16) 

Hence by {^J^ and (|1|) 

-A(v^ + s\x\-\ij - ^) < 

for |x| > ri. We fix M such that 

?/'(ri) < ^{ri). 

li ip ^ for |x| — ^ oo, we immediately get 

ip <ip for |x| > ri (2-17) 

from the maximum principle. The fact that J ip"^ < oo and ip G L"^^ contradicts 
our assumption. Unfortunately, we only know that ip ^ oo as \x\ — i>ooina weak 
sense, namely ■?/' G L^, so the authors in [BBL] used a variant of Stampaccia's 
method to verify the statement of the Lemma, which also works in our case. □ 



Mimicking the proof of |[BBH| Lemma 10 and using the fact that t'^{iI)'^)%P — (pip 



is continuous but not differentiable we get 

2.7. LEMMA. ip > {) everywhere and ip & C"^ , except at x = 0. 

Using ( p^.2D in the proof of | |BBL| ] Lemma 11 and afterwards following the 
proof of Lemma 13 we additionally get 

2.8. PROPOSITION. A^^ = Jip^ > Z. 



Before concluding the proof of Theorem 2.1, we need a final lemma, which is 



the equivalent to |PBL|| Lemma 14. 



2.9. LEMMA. For every N > we have inf{£^™[p]|p e D and J p = N} 
mi{£^^^'^[p]\p eDandJp< N}. 



Proof of Theorem 
For every we set 



E{N) = inf{^^™[p]|p G Dand J p< N}. 
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Obviously E{N) is non-increasing and convex. The same proof as in Proposition 
|2.3| shows that there exists a pn & D with J Pn ^ N and 



Cj ; 



g^''''^^[p^] = E{N). 

With Nc = Jtp'^ it is clear that E{N) is constant for > A^^: E{N) = E{N, 
while E{N) is strictly decreasing on the interval [0, Nc]. For N < : J p^ = N, 
which implies that ( |2.^ ) has a unique solution. On the other hand we deduce 
from Lemma p]0 that for N > (|2.3|) has no solution, which concludes the proof 



of Theorem 2.1. □ 



After having guaranteed the existence of a minimizing density p at for (|2.3| ) , we 
can derive an Euler-Lagrange equation under the variational restriction J p = N. 

1/2 

2.10. PROPOSITION. Denote ip = pj^ , with p^ the minimizing density for 



inf{£™[p]|pGA/p = iV} 



under the restriction N < Nc- Then we have 

- + r^(^2)^ -v>^ = p{N)^, (2.18) 

where 

KN) = ^E{N). (2.19) 

Proof. The derivation of ( p^.l8| ) works analogously to ( p^.ll| ) apart from the dif- 
ference that p is the Lagrange parameter for the restriction J p = N. We can 
infer 

"^"^'^ KmJipM-p), (2.20) 



^f™[tp + (l-t)p;v] 



dt 

which implies by means of convexity of the functional 

^™[p^] - f ^™[p] > p{N) I (p^ - p), 

or equivalently for every A^' 

E{N) - E{N') > p{N){N - N'). (2.21) 
On the other hand we derive from ( p.20|) 

£:™[tp + (1 - t)p^] - ^^™[p^] = t/i(A^) I (p^ - p) + o(t), 

which yields with p = 2pAr and p = jPn, respectively, 

E{N ± tN) - E{N) < ±p{N)tN + o{t). (2.22) 
Hence (g]^ and (p:^ ) together imply (gJ^. □ 
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Next we take a look at the behavior at infinity of the minimizing densities pn- 
At least for N < Nc one gets exponential decay. 

2.11. PROPOSITION, (a) Let n< 0, which is equivalent to N < N^, then 
for every 6 > 0, with fi < —6, there exists a constant M, such that, for the 
corresponding minimizer ip = p]^"^ , 

iP < Me-^''"l"l. (2.23) 

(b) Let N = Nc, then for every 6 < Nc — Z there is a constant M, such that 

^ < Me~^^'\^\^''\ (2.24) 

Proof. Note, that we have not yet shown that ip ^ as \x\ ^ooina strong 
sense, for we only know ip & L"^ . From equation ( p.l8| ) we derive —Aip < Vi/), 
which implies (—A + /)■?/' < {V + Since we know (V^ + /)■?/' G L^, recall that 
{V + I) e L'^ + L°° and ?/; e fl we conclude from 

^<(_A + /)-i[(l^ + /)^], (2.25) 



that Tp —y Q infinity, e.g. from the well known fact ( ||LS|| Lemma 11.25) that 
the convolution f * g oi two functions / G L'^'^g G L^, with 1/p + 1/g = 1, goes 
to at infinity. 

(a) (cf. [LI] Theorem 7.24.) Let 5 < -p. From (|218|) we get 

{-^^5)ij=[-T'B{ilj^) + ^ + p + S\ij (2.26) 

and 

7/- = (-A + 5)-\-t'b{^^) + V + P + 5]^. (2.27) 

Since V,ip ^ Q as \x\ — * oo, there is a ri, such that [— t^('?/'^) + + p + 5] <0 
for \x\ > Ti- This implies 

^Pix)< I {A7t\x-y\)-^e-''^'\^-y\{[-T'B{i,^) + ^ + p + 5]i,){y)dy < (2.28) 

J\y\<ri 

and ( p. 231) with 

M = supe'^'''^'^ / {An\x-y\)-\{-r'j,{,lj^) + ip + p + 5]ij){y)dy, (2.29) 

X J\y\<ri 

(b) This follows directly from ( p.l7| ) and the fact that t/^— s>Oas|x|— >oo. □ 

We state a final proposition concerning the behavior of the chemical potential 
at = (which is — oo in the usual TF theory), because of the simple and 
illuminating proof. 
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2.12. PROPOSITION. Leteo = ~\Z^ he the smallest eigenvalue of the Schrd- 
dinger operator —A — Z/\x\. Then 

'' = e„. (2.30) 



N=0 



Proof, (cf. [BBL] Lemma 15.) Let <f{x) be the normalized eigenvector of —A — 
Z/\x\ belonging to the lowest eigenvalue Cq = —\Z'^ and let pat = NLp{xY. Then 

E{N) < £:^™M = j \Vpf\'' - jypN + j rB{pN) + D{pN,pj,) 

= A^[(v9, -Av9) - Z(v9, |x|"V)] + j TnipN) + DipN,pN) 
< Nco + CiN^/^ + C2N\ 
On the other hand we have 

E{N) > inf^S^j \Vp^/^\^ - j Vp 



which altogether implies 



A^inf spec{-A - Z\x\'^} = Ne. 



y E{N) 

lim -— = Cq. 

N 



Taking into account that E{0) = this is equivalent to ( p.30|) . □ 



2.2 The Scott correction 

If one takes a look at Lieb's proof [LI] that in the usual TF theory without mag- 
netic fields the von Weizsacker term produces the Scott correction, one realizes 
that the main correction comes from distances of order Z~^ from the nucleus. 
It is thus reasonable to guess that in the MTF theory the von Weizsacker term 
produces the Scott correction as long as p^^^, the density corresponding to the 
MTF energy ( |1.6| ), is well approximated by the usual TF density p"^^ up to dis- 
tances of order Z~^ from the nucleus. This condition is equivalent to the demand 
that t'^{p) is proportional to p^/^ for r ~ Z~^ and this is the case for B <^ Z^. 
In other words, for B <^ Z^ the von Weizsacker term produces a Z^ correction at 
the distance of order Z~^ from the nucleus. At the edge of the MTF atom, the 
radius is known |[LSY2|| to be proportional to Z^I^B'"^!^ and the lowest Landau 



band is occupied, which leads to r'^^p^^^^ = 4:7i^B~'^{p^^^Y in the outer region. 
Hence, one computes very easily, by using p^"^^ as comparison density, that the 
correction coming from the edge of the atom is of order B^^^Z^^^. 

Thus, for B > Z"^^^ the correction from the edge of the atom overcomes the 
Scott correction and Theorem |1.4| follows by using the variational density ( |2.35| ). 
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So it suffices to prove Theorem |1.4| for B < Z"^!^. 
Proof of Theorem [IT 



First of all notice that there is an rs, such that for r > tb the density p^'^^ 
corresponds to the lowest Landau band, i.e. 

r^(p^^"W) = ^(P^^^(r))^ for r > r^, (2.31) 

and for r < rs we have 

n.ip^'^'ir))'/' < rUp^^^(r)) < n,ip^^^ {r))'/^ (2.32) 



with K3 = 0.83(37r^)^/^. Using ( p^.31 ) and ( |2.32| ) one realizes that if one fixes any 



e>0 with B < Z^'" there is a 5 > such that > Z^'^. 
Lower bound: 

We know from Section 2.1 that there exists a po satisfying 

^MTFW ^ ^MTFW[^^] ^ J + J ^^^1/2^2 _ J y p, + D{p„p,). (2.33) 

Denote Z\x\-^ = V = V + H, with 

H = Z /r ~ Z"^ /h for r <h/Z and otherwise, 
= Z'^ /h for r < h/Z and Z/r otherwise. 

Now let us rewrite the energy functional in the following way: 

^™[po] = / Mpo) - J pT^pT - I pTVpT -Jhpo + D{po,po) 

Observe that -A - H > Mp{{V p^^^ ,V p^^'^) - {p^'^ ,Hp^/^)}, which by using 
Sobolev's inequality can be bounded from below by 

-A-if>inf{||p||3-||p||3||if||3/2}. 
p 

Since || H ||3/2~ h we can guarantee —A — H > with b small enough. Choosing 
such a 6 we derive 

^MTFW > g^^^[po,V] 

> E''^^[V]= [ Tip) - [ pV + Diprp) 
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The density p minimizes S^'^^lp, V] and therefore fulfills the TF equation: 

^b(p) = V - p*\x\'\ 
By means of (1021) we get that J Hp = 0{Z^) for B < Z'^l^ which yields 

^MTFW>^MTF^Q(^2)^ (2.34) 

Upper bound: 

In order to get an upper bound we use a variational density p in the following 
way: 

r f^iZ-') for r < Z-\ 
p(r) = i pTF(^^) < r < rs, (2.35) 

where p"^^ indicates that we essentially use the usual TF density p^^ . Only in 
the region eZ/B < r <rB we eventually have to modify p"""^, such that 

lE^T^F -£^^^[p]\<OiZ^). 

(This e.g. can be done by following the way of and taking the modification 
W of the effective potential 0^"^^ and defining ATcp^^ = A{W - Z\x\~'^).) By 
means of our variational density p we get 

UjVp'/r = OiZ'), (2.36) 
Ir>r, |VP^/T = 0{B'/'zy% (2.37) 

For B < Z'^l^ this leads to 

^MTFW<^MTF^Q(^2)^ 

which together with ( |2.34| ) completes the proof of Theorem |1.4| . □ 



3 Gradient corrections for STF type theories 

3.1 The functional ([T^^) 

Starting point in this section is the functional 

^'^[p] = ^,lp' + j-,li'^p'^r-lvp + Dip,p), (3.1) 
with V and D{p,p) defined as in (|1.3|). Compared to ( |1.50D we rewrite 
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and for simplicity forget about the numerical constant, which does not effect any 
mathematical statements. The corresponding energy is defined as 

E{N,Z,B) = M{S'^[p]\p eD,jp<N}, (3.2) 
where the domain D is given by 

D = {p\p > 0, p e n L^ ^ 

In analogy to Section 2 we first consider the problem 

Mm{E'^[p]\pED}, (3.3) 

with 

D = {p\p >0,pe L^ G L^, D{p, p) < oo} (3.4) 

and show that the minimum is achieved for a unique po. By means of the cor- 
responding Euler-Lagrange equation we shall deduce that po is in L^(]R^), more 
precisely / po = Z. 

First of all, we collect some properties of (|3.1| ). 

3.1. LEMMA. There are positive constants a, C , so that 

^'^[P]>«(llp||3 + j rB(p) + ||V//2||^ + D(p,p))-C, (3.5) 



Proof. This is a consequence of Lemma 2 in | BBL |, which tells us that for every 
e > there exists a constant Ce so that 



Vp<e\\ph+C,D{p^pfl^ 
for every p > 0. □ 

3.2. LEMMA. £^'^[p] is strictly convex m p. 

Proof. This follows immediately from the strict convexity of (Vp^^^)^ and 8^"^^ . 

□ 

By means of these Lemmas we can prove the existence of a unique minimizer in 
D. 

3.3. PROPOSITION. The Minimum of £^^[p] is achieved by a unique po E D. 
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Proof. The proof is similar to that of Proposition 2.3. Let p„ be a minimizing 



sequence. By Lemma we have 



||Pn||^<C^, ||Vpf||^<C, D(p.,p„)<a (3.6) 

With Banach-Alaoglu theorem we therefore extract a subsequence, still de- 
noted by pn, such that 

Pn Po weakly in L^, (3.7) 
Vpn^^ ^ Vpo^^ weakly in L^. (3.8) 

Furthermore pf/'^ is bounded in if^, which implies that there exists a further 
subsequence, again denoted as pn, with 

3/2 3/2 

Pn Po a. e.. 

(This relies on the fact that for a smooth bounded domain Q, H^{Q) is relatively 
compact in L^{Q).) Hence, using Fatou's Lemma we get 

liminf D(p„,p„) > D{po,po), 

and by the weak lower semicontinuity of L^-norms we deduce 

liminf/(Vpr)^ >/(Vpf)^ 
liminf /p3 > Jp3 

In order to prove J Vpn — > / Vpo, we decompose V = V1 + V2 such that both 
functions are in C°°. With Vi G L^/^, { ^7f\) implies 



VlPn -* / ViPq. 



On the other hand V2 fulfills 

J V2Pn = J l^2[-A(p„ * \X\'')] = J {-AV2){pn 

which converges to 



* \x 



-AV2)ipo*\x\-') = J V2P0, 

for — AV2 G L^/^ and || pn * He is bounded. 
Thus 

limmf£'^[pr,]>£'^[po]. 

The uniqueness is an immediate consequence of the strict convexity of the func- 
tional. □ 
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For this minimizing po we can derive an Euler-Lagrange equation. 
3.4. PROPOSITION. The minimizing V^^/s = p^/^ satisfies 

{-A + B)^ = ^^r''\ (3.9) 

with ip = V ~ ip"^^"^ * m the sense of distributions, on the set where if) > 

Proof. The uniqueness of the minimum, the spherical symmetry of the functional 
( |3.1|) and the fact that ip G implies the continuity of ip away from the origin. 
With y9 G L^Q^ the Equation (|3.9| ) has a meaning in the sense of distributions on 
the domain {x\ip{x) > 0}. Consider the set 

D = {r]\r] G H\ D{rf'^ ,rfl^) < oo}. 

If ?7 G Z), then p = {ri^y/^ G D and 

For all ?7 G i) we find (t){ip) < (t>{f])- 

Let i G C^, then using + ^Olt=o = we infer 

□ 



3.5. PROPOSITION. is hounded andip is in C°° away from the origin and 
eventual points with ip{x) =0. 

Proof. Denote 

VL, = {x\ip{x) > e}. 

On this domain we have 

(-A + 5)^ = /, 

with / G Lf^^, since (p G L^^^. So we conclude from standard elliptic arguments 
(e.g. [LL] Section 10) that is bounded, hence continuous everywhere. From 

= 47r(V'^/^(x) - Z6{x)) 

we get the two times differentiability of ip away from the origin and as long as 
ip > 0. By means of ( p.9|) and a standard bootstrap argument we conclude 
G □ 



H November 29, 2000 



24 



3.6. PROPOSITION. / V^^/s = z. 

Proof. Suppose by contradiction / ip"^^^ = \ ^ Z. We do not assume A to be 
finite. Then, as in ( |2.14| ), we get that there is some r\ and an e > 0, such that 

|</.(x)| >e/r, (3.10) 

for |x| > r\. Since il) G L^, we conclude < M/r^/"^ for r large enough. Hence, 
from ( p.9|) we get 

|(-A + 5)7/^1 = > er-^/^ 

which is a contradiction to G H^. □ 



3.7. Remark. The Proposition p.6| is interesting, since gradient corrections usu- 
ally give rise to binding of additional electrons. The reason that this is not the 
case in relies on the negative potent of ip on the right side of ( |3.9| ). This fact 
forces the potential (p to fall off much faster than 0(l/r). 

Since we only consider atoms with a point nucleus, we get the following remark 
for the minimizing po- 

3.8. PROPOSITION, po is a symmetric non increasing function of |x|. 

Proof. As we have already argued above, the symmetry of po follows from the 
symmetry of the functional and the uniqueness of po. 

Denote the non increasing rearrangement of po- (For definition see e.g. 
[[LL|| Section 3.3.) From the fact that j po < Z and [LI] Theorem 2.12 we get 
^^'^^[pl] < S^'^^Ipo]- iS Lemma 7.17 implies 



|V(pf)*p> y |V(pf)p (3.11) 
and again from |LL| 3.3 (v) we get {pl^^y = (p*)3/2, which proves the statement. 



□ 



3.9. PROPOSITION. ip has compact support. 

Proof. Inserting Aip = into (|3.9| ) yields the following equation for the po- 
tential, away from the origin: 

(A^)V2[A(Av.)^/2] = (A^)2-^, (3.12) 

where we replaced the constants by one. Since is spherical symmetric we can use 
the ansatz (y9 = \x{r) and obtain by (p.l2|) the following fourth order equation: 
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= I[^X' - ^XX!" + {Ix")'] (3.13) 

We can see that in the surrounding of each point tq > there exists a local 
solution 

X = aeir - rof + 0((r - roY), 

with flg = (const. )rQ. Since x = is also a solution of ( |3.13|) , every composed 
function 

X = cbeir — Tq)^ + 0{{r — roY) for r < Tq and X = for r > Tq 

is a local solution around tq. 

Away from the solutions can be uniquely continued up to r = 0. Hence, 
there exists a solution x and a ri > 0, such that x(0) = Z and suppx = [OjTi]. 

Repeating the argument of [LI] Theorem 2.6, one can show that a solution of 
(13:91), with e and / V'^/^ < oo, (ip'^Y^^ uniquely determines the minimum 
of the functional ( p.l| ). Therefore x{f) uniquely determines the selfconsistent 
potential (f) = x/f, which implies that ijj = pl^"^ = {A(j)Y^'^ has compact support, 
too. □ 



3.10. Remark. The preceding three propositions are equivalent to those for the 
STF theory. This confirms that the p(Vp)^ term only amounts to changes close 
to the nucleus. 

By the convexity of the functional (|3.1|) one easily derives the following properties 
for the energy E^'^{N, Z, B): 

3.11. PROPOSITION. E^^{N, Z, B) is convex as a function of N and strictly 
monoton decreasing on the interval [0,Z]. For N > Z we get E^'^{N,Z,B) = 
E'^iZ,Z,B). 

Next we prove Theorem p..7| . 
Proof of Theorem \1 . 1\ : Upper bound: 
We use the comparison density p, with 

= / P'^^(^i^) for r < Is = B-'/^ 
1^ p^^^(r) otherwise 

and immediately get 
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Lower bound: 

Let p be the minimizer of the energy (|3.2| ), for given B, Z, N. We can rewrite 
" ''p' + ^Jl'^p'^'l'' [hp- fvp + D{p,p), 



J B\ 

where H = ^ — ^-j-^ for r < Mb and otherwise, V = ^-j^ for r < Mb and 
Z/r otherwise. Looking at the term 

ij / IVp^'T - / Hp, 
we infer by means of the Sobolev inequahty the estimate 

i_ f |Vp3/2|2 _ /"^p > _i_ II p ||3 _ II ^ ll^^^ll (315) 



Since || H ||9/8~ b^^^^ we can choose b such that (|3.15|) > 0. Hence 



K 



E'^ > p'- J Vp + D{p,p) 

> E^^^[V]=£^^^[V,p] 



Since / pH = 0{Z^^'^B^^^) we prove the proposition. □ 

We see that the main contribution of the correction (|L50|) comes from the 
radius Ib = B~^^'^, in other words the gradient correction repairs the infinity 
of the STF density at a distance from the nucleus. This infinity stems 

from the fact that in the STF theory the full potential is involved, although 
particles in the lowest Landau band, which are smeared over a radius of at least 
never see the full Coulomb singularity. 

Moreover, the Tomishima-Shinjo correction (|L46| ) orthogonal to the magnetic 
field remedies the singularity of the Coulomb potential in a similar way as the 
isotropic gradient term. We will see that the same effect, as caused by and 
(|1.46D, is also naturally obtained by using the DSTF functional. 



3.2 A discrete von Weizsacker functional 
3.2.1 The DSTF functional 

First of all we are going to collect some information about the DSTF functional 
(which are rigorously proved in a companion work). The DSTF functional 

^DSTF^^] =J2Uf pU^) - I V^{z)p^{z)dz\ +Dip,p), (3.16) 
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with Vn and Vm,n as in ( |1.35| ), is defined on the domain 

D = {p\ J2 J Prn < 0O,J2 J Prn < ^, D{p, p) < Oo} (3.17) 

m m 

with corresponding energy 



p G D and ^ /"p^ < AT I . (3. 

m ) 



Following the considerations of ||LSY2|| one can easily see that ( p. 161 ) is convex 



and bounded from below on = {p| p G -D, I Pm ^ N} and derive the 
following Theorem: 

3.12. THEOREM. With N < Z fixed there exists a unique minimizer for 
£^^^^, under the restriction Em I Pm < N, i.e. E^^^^{N,Z,B) = S^^'^^lp^]. 
Moreover, p^ satisfies / p^ = N. 

Furthermore each minimizer p^ obeys the coupled TF equations 

3K{pZ{z)f = [ZVUz)-J2 I Kn,„(^-^)P^(^)+M^^)]+ V(mGNo), (3.19) 

where p{N) is the Lagrange parameter belonging to the restriction J pm = N 
and [] + , with = t for t > and = otherwise, corresponds to the fact 
that the functional is only varied over positive functions, i.e. Pn > Vn G Nq. 
By means of the notation 

ZVUz) ^-."(^ - ^')Pn{^') + f^{N) = ^i?(^) (3.20) 

and inserting in the TF equation we can rewrite the energy E^^"^^ {N, Z, B) = 

E^'^^iN,Z,B) = Y: I l[p'-V^\z)]-^^ + Np-b(p^,p-). (3.21) 



3.2.2 Semiclassical approximation of £^conf 

First of all we want to state a useful theorem concerning the sum of the negative 
eigenvalues of the one-particle operator HqIHa + V^jHo, where Ha = ((— + 
A(x)) ■ 0")^ and (p is an axialy symmetric potential ip{r, z) with r = |x_l|. 
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3.13. THEOREM. Let (f — (p{r, z) be axially symmetric. Then one can write 
the trace of the negative part of the operator IIo[Ha + f\Ilo as the sum of one- 
dimensional traces, i.e. 



with 



meNo 



(3.22) 



(3.23) 



Proof. Let denote the angular momentum operator parallel to the magnetic 
field. Since (f — (p{r, z), we have 



(3.24) 



which implies that the eigenvectors of the operator no[i?A + <^]no are of the form 
\m,i) = $^(x^)/4(^). 

Hence we can write the sum of the negative eigenvalues as 

Tt[Uo[Ha + <p]UoU = Y.{m,i\[Uo[HA + ^]Ilo]\m,i) 

m,i 

= EfE(/-[-^' + <^^"^n*j), (3.25) 

showing that the /^'s are the eigenvectors of the one-dimensional operator — 9^ + 

cp^'^^z). 

3.14. PROPOSITION. Let N, Z, B be fixed. Then 

e2,AN,Z,B) < E''^^^{N,Z,B) + R, + R2 + R3, 
E?oAN:Z,B) > E^^^^{N,Z,B)-R,-cJpt/\ 



□ 

(3.26) 
(3.27) 



with 



Ri 

R2 

R3 



J2 

meNo ^ 

XN<>'i<K^) 




dpdz^ 



(3.28) 

(3.29) 
(3.30) 



where Aj and p are defined in the proof. 
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Proof. Upper bound: 

First of all we note that for any comparison wave function and fixed integer 
we get 



N 



(3.31) 



1=1 



with 



If we set p(x) = J2m\^rn{'^±)\'^Pmi^)^ subtract p * |x| ^ — p{N) in 



the scalar product and use ip 



m 



A ... A (f)N, where 0i is the eigenvector 



corresponding to the i-th lowest eigenvalue A, of the one-particle operator 

UoiHA - Z\x\-' + p * |x|-i - /i(iV))no, (3.32) 
as comparison wave function, ( p.31|) reads 



E 



conf 



< Ti[UoiHA - Z\^\-' + p * |x|-^ - p(iV))no]- - 2D{p^, p) 
+D{p^,p^) + Np{N)+ J2 (^^-MA^))- 

A]v<Ai<At(Af) 



(3.33) 



Applying Theorem 3.13 and (|3.21 ) to the inequality above, we finally arrive at 
the upper bound ( |3.31D . 
Lower bound: 



Let tjj denote the minimizer of (|1.22| ), i.e. tjj = V'conf- So after again adding and 
subtracting p * |x|~^ — p(A^) and using the Lieb-Oxford inequality ||L0|| , we can 
write the lower bound on -E^^f as follows: 

^clf = {^,KhnK) 

N 
i=l 

+D{p^, p^) + NpiN) - 2D(p^, ~P)-cJ pf 
> Tt[Uo{Ha - Z\^\-' + p * Ixj-i - p{N))Uo]- + Np 

-D{p,p)-C [ pf (3.34) 



Using ( p.21| ) we arrive at ( |3.27| ). 



□ 
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3.15. Remark. Due to ( |3.26| ) and ( |3.27| ) the main contribution to the difference 
between E^,^^ and E^^'^^ is given by 



Ri 



dpdz 
2n 



(3.35) 



which, by definition, shows that _£;dstf jg ^^^q natural semiclassical approximation 
of i^lf ■ 



As a corollary, for example by following the way of ||LSY2|| and using coherent 
states and the above estimates, or simply by estimating the difference between 
E^^'^^ and E^'^^ one gets 

3.16. COROLLARY. If Z ^ oo with N/Z fixed and B/Z^ 0, then 



E^^^,{N, Z, B)/E^^'' (AT, Z, B) ^ 1. 

The DSTF theory is equivalent to a three dimensional functional using mod- 
ified Coulomb potentials, 

K,(x)= J]x"(xx)KW (3.36) 

n 

replacing the attractive Coulomb potential and 

Vn,m{z - z')Xn(x±)Xm(x'_L) (3.37) 

replacing the Coulomb repulsion, with 



1 for ^f{2n)/B <\^i_\< ^2{n + l)/B , . 



otherwise. 

I.e., for the respective minimizing density of a resulting MSTF functional, we have 
^MSTF(^) _ M. '£^p^'^^{z)xmM as Well as E^^tf _ ^dstf ^^le energy. 
Since V^{Q) ~ S^/^, 1/^(x) can be regarded clS cl cut off Coulomb potential 

= I fo,. > (3.39) 

Hence, it is obvious that the main contribution to the difference E^^"^^ — E^^'^ 
stems from the Coulomb potential in the region r < given by the term 

B I |(^STF|3/2^ 0(^1/4^3/2^^ (3_4Q) 

io<r<B~V2 

which leads to the relation 

E°STF(jv, Z, B) - E^^^'iN, Z, B) = 0{B^'^Z^'^). (3.41) 

The comparison with Theorem [l.T] shows that the DSTF theory has the same 
effect as the introduction of the gradient correction in ( |1.5CI| ) as well as in ( |1.46| ). 
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3.3 The discrete von Weizsacker functional 

The variable of the DSTF functional is given by a sum of one dimensional densi- 
ties, emphasizing the character of lowest Landau band particles, whose positions 
orthogonal to the magnetic field are "frozen" and they therefore only move par- 
allel to the magnetic field. Now taking into account the result (cf. e.g. Shao 
that the first order correction to the semiclassical description of the one- 



dimensional free Fermi gas is given by —(1/3) J (dz\/p{z))'^, we are motivated to 
propose the already mentioned discrete von Weizsacker functional 

S^'^IP] = E (-^/ + ^ J pUz)-Z j V,^{z)pm{z)^ +Dip,p). 

(3.42) 

Since the von Weizsacker term appears with negative sign, ( p. 42 ) has the same de- 
fects as the Tomishima-Shinjo functional ( |1.46| ), i.e. it is not bounded from below 
and not convex. Precisely these two features (boundedness and convexity) of the 
semiclassical TP functionals, as well the (M)TPW functional, provided not only 
the existence of a minimizer but the existence of a solution of the corresponding 
TP equation. 

As a way out of this problem we can define the energy corresponding to ( |3.42| ) 
by means of stationary solutions p^ , whose variational derivative vanish under 
the restriction J p^ = N and p > 0. In order to avoid the assumption 
of positivity we concentrate on real functions ip, with ip"^ = p and consider the 
functional 



m = E (-^/ \d^i^m{z)\'' + K j i^Uz)-Z j Vrn{z)i^l{z)^ +D{,p^r). 

(3.43) 



mSNo 

Let 



^ = {^1 E < ^' E < ^' E / (^^^-)' < (3.44) 

m m 

be the domain of ( p.43| ). The question for stationary points in D, under the 
restriction J ip"^ = N, is equivalent to the existence of a Lagrange parameter 
p{N) and a ip'^ , so that 



I { £[i;^ + tr^] + piN) li^P'^ + tvY 



0, (3.45) 



i=0 

for each rj E D. (|3.45| ) yields the Euler-Lagrange equation, denoting ip^ = ip, 

{l/3)dl^Uz) + 3K^i(z) = [<^y (^) - p{N)]tPUz) Vm G Nq. (3.46) 

Starting from ( p.43| ), (|3.46|) a priori only exists in the sense of distributions, but if 
there is a solution for (|3.43|) then one can conclude that it is even in C°°(]R\ {0}). 
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If there is a solution ip^ for G [0, Nc], with Nc > Z, then we can define the 
corresponding energy by 

E°W(Ar, Z, B) = £[ilj^] = (3.47) 

3.3.1 Recovering the exchange term 

Now we even go a step further. Following the reflections of Section 1.3.1 con- 
cerning the magnitude of the negative von Weizsacker term we guess that 8^"^ 
is equivalent to the DDM functional ( |1.33| ) for B <^ Z^. Hence looking at ( |1.42| ) 
we suggest another functional, 

= E„.eMo (-(1/3) / \d.^l^? + ^IpL{z) 

- ZJVUz)pU^) - cHB / Z^l^) j pI) + 5(p,p), (3.48) 

where we recover the exchange energy, which could be compared with the Thomas- 
Fermi-Dirac-von Weizsacker functional ([LI]) in the 5 = case. 

3.3.2 The one-dimensional DW functional 



In this section we study a toy model obtained by reducing (|1.52|) to a one 



dimensional functional and dropping the Coulomb repulsion, which leads to the 
functional 

£i°W[p] = -\j |9,v^P + ^ j P'{z) -Z j V,{z)p{z). (3.49) 
First of all we consider the corresponding TF functional 

£Z\P\ = ^ [ P'(^) -Z [ Voiz)piz), (3.50) 



for which we easily get the following lemma: 

3.17. LEMMA. £^]j°^[p], with p e L^(M.), is bounded from below, and there 
exists a unique minimizing density po G L^, with SjP^lpo] = E^^{Z,B). 

Proof. Since J Vop <\\ Vq II3/2II p Us, we get 

4°Z[P]>IIP|I3-|IK.||3/2||P||3. (3.51) 

Minimizing over || p II3, we get the first part of the lemma. The proof of second 
part works analogously to Propositions 3.2 and 3.3. □ 
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For the minimizing density po we get the simple TF equation 

SkpKz) = ZVoiz). (3.52) 
By the definition ( p,.35| ) one easily sees the relation 

Vo{z) ^ = B'/X\B'/h), (3.53) 
which implies an interesting scaling relation for the energy E^^{Z, B): 

3.18. LEMMA. 

E'^{Z, B) = Z'^/^B^/^E^^il, 1). (3.54) 
Proof. Using the scaling relation ( p.53| ) and defining 

p{z) = B^'^Z^'^p{B^l^z), (3.55) 

we get 

£Z\P] = B"'Z''^£I^M- (3.56) 

□ 

In the next Theorem we point out that for B Z^, the energy ( |3.54|) is 

the semiclassical approximation of Tr[— 9^ — ZVo{z)]^, the sum of all negative 
eigenvalues of —d"^ — ZVo{z). 

3.19. THEOREM. Let B < Z'^ and ip e C^{B{Q, B-^/'^)). Th 



en 



Tr{^p[-dl - ZVoiz)].) = Z^'^B'I^ j j ^i.[p^ - V^^z)]. - O^B^'^Z^'^), 

(3.57) 

and there is a constant C, such that 

\TT[-dl - ZVoiz)]- - Z'^'^B^'^E^^il, 1)1 < CB'^'^Z^'^ (3.58) 

Proof. Let us rewrite 

- dl - ZV^^iz) = B^'^Z[-iB^'^Z)-^dl - V^iB^/'^z)] (3.59) 

and define the unitary operator 

iUil)i))iz)=l^/^i)ilz). (3.60) 

With / = B^^/"^ and the fact that unitary transformations do not change traces 
we derive 

TT[-dl - ZVo^iz)]_ = B'/'ZTT[-h'dl - V^iz)]_, (3.61) 

where 

h = B^'^Z-^/^. (3.62) 
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So we have the self-adjoint Schrodinger operator 

H = -h'dl-V,\z), 

and its symbol 



Using the relation 



h{z,p)=p'-V,\z). 



d^Vf^'iz)] < C on 5(0,2) = {z\ \z\ < 2} 



(3.63) 
(3.64) 

(3.65) 



and by means of some cut off function ip G 



Theorem 6.1 yields 



TT{^{z)[-h'dl-V,\z)]^) = h 



-1 




dpdz 
2n 



tP{z)[p'-V,\z)]^-0{h), (3.66) 



which together with (|3.61| ) and (|3.62|) immediately implies (|3.57|) . 
Recall, that by = t for t < and otherwise. 

In order to tackle the outer zone {z\l < \z\ < oo}, we note that the potential 



fulfills 



Id^Vo'iz)] < Clz] 



for \z\ > 1. 



(3.67) 



So by definition of the scaling functions l{z) = \z\ and f{z) = 1, we can use ||^] 
Theorem 7.1, which states that with a cut off function ip G C°°{B{1, oc)), we get 

dpdz 



TT{^[-h'dl-V,\z)]^)-h-' 




27r 



-^P[p'-V,\z)l 



< h I dzl~^ < h. 

(3.68) 

Combining (|3.66| ) and (|3.68| ) together with ( |3.61 ) and ( p. 62 ) completes the proof 
of Theorem |3.19| . □ 

Turning back to the functional ( |3.49| ), we recall that the corresponding energy 
has to be defined by means of the solution of the TF equation 



dl^{z) + 3k^^{z) - ZV^{z)^{z) = 0. 
Introducing the scaling relations 

ij{z) = B'I^Z'l^ip{zB'l''), z ^ zB^/\ 
(|3.69|) can be written as 

^dl^{z)+?,K^\z)-V,\z)=Q. 



(3.69) 
(3.70) 
(3.71) 



We set e = B^l"^ jZ -C 1 and make the ansatz ip = Lpo + eipi to get an approximate 
solution of ( p.71| ). Rescaling and inserting into ( p.49|) yields 



E 



IDW 



(3.72) 



which is in accordance with (|3.57| ) and justifies a posteriori the introduction of 
the negative von Weizsacker term in ( p..52| ) and ( |3.49| ). 
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4 Some concluding remarks 

Throughout this paper we have studied, among others, two functionals, the 
MTFW, which represents an approximation to the full quantum mechanical en- 
ergy E^, and the DW functional, which should approximate the ground state 
energy -E^^f of particles in the lowest Landau band. Now we can ask for the 
magnitude of B, for which the use of DW becomes more reasonable than MTFW. 
This question is connected with the estimation of 

E<iiN,Z,B)-E2,^,{N,Z,B)\. (4.1) 

For B ^ this can be compared with the difference of the corresponding 
semiclassical approximations 

B J |0^'iTF(a;)|3/2^a; + 2B / - 2tBffdx -Bj {x)ff dx 

<CJ \Z\x\-^ - 2B\^^dx < CZyB^/\ (4.2) 

From the preceding sections we guess, on the one hand, 

lE"""^ -E2,J<o{B^/'^Z'/'), (4.3) 

and on the other hand we know 

I^MTFW _ < Q^z^^ ^ 0(5^/5^3/5). (4 4) 

Hence, one might expect the DW theory to give the better description of the 
quantum mechanical energy, if Z^/B^^"^ < B^I^Z'^I^^ or, in other words, B > 

^24/13 _ 

Summing up, we suggest to use MTFW theory for B < Z^^/^^ and DW for 

^24/13 <B^Z^ 
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